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Abstract 

To clarify the role of potassium inwardly-rectifying-channel, subfamily-J, member 11 {KCNJll) variation in susceptibility to 
type 2 diabetes {T2D), we performed a systematic meta-analysis to investigate the association between the KCNJll E23K 
polymorphism (rs5219) and the T2D in different genetic models. Databases including PubMed, Medline, ElVlBASE, and ISI 
Web of Science were searched to identify relevant studies. A total of 48 published studies involving 56,349 T2D cases, 
81,800 controls, and 483 family trios were included in this meta-analysis. Overall, the E23K polymorphism was significantly 
associated with increased T2D risk with per-allele odds ratio (OR) of 1.12 (95% CI: 1.09-1.16; P<10"^). The summary OR for 
T2D was 1.09 (95% CI: 1.03-1.14; P<10"^), and 1.26 (95% CI: 1.17-1.35; P<10"\ for heterozygous and homozygous, 
respectively. Similar results were also detected under dominant and recessive genetic models. When stratified by ethnicity, 
significantly increased risks were found for the polymorphism in Caucasians and East Asians. However, no such associations 
were detected among Indian and other ethnic populations. Significant associations were also observed in the stratified 
analyses according to different mean BMI of cases and sample size. Although significant between study heterogeneity was 
identified, meta-regression analysis suggested that the BMI of controls significantly correlated with the magnitude of the 
genetic effect. The current meta-analysis demonstrated that a modest but statistically significant effect of the 23K allele of 
rs5219 polymorphism in susceptibility to T2D. But the contribution of its genetic variants to the epidemic of T2D in Indian 
and other ethnic populations appears to be relatively low. 
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Introduction 

Type 2 diabetes (T2D) is a complex metabolic disease resulting 
from reduced insulin secretion and peripheral insulin resistance. 
By coupling cell metabolism with membrane potential, adenosine 
triphosphate-sensitive potassium channel (KATP) play a central 
role in regulation of insulin secretion in pancreatic-P cells. [1]. The 
KATP channel is a hetero-octamer of K"*" inward rectifier Kir6.2 
[KCNJll) and regulatory sulfonylurea receptor SURl subunits 
[aBCCS] [2]. Mutations in both KCNJll and ABCCS cause 
neonatal diabetes and congenital hyper-insulinemia in humans 
[3,4]. In addition, KCNJll gene knock-out mice are characterized 
by defects in insulin secretion in response to either glucose or 
tolbutamide [5]. 

As a candidate gene for T2D in humans, a nonsynonymous 
E23K variant (rs5219) which results from a G — > A transition in 
codon 23 in the NH2-terminal tail of Kir6.2 was identified [6]. 
With spectacular advance in genotyping method in recent years, 
larger-scale genetic association study concerning the relationship 
between the E23K polymorphism and T2D susceptibility has been 
conducted in various populations. However, inconsistent results 
have appeared in the literature. Such inconsistence may be due to 
chance, insufficient power of limited sample size, or bias in study 
design (e.g., inappropriate control selection). Alternatively, these 



disparate findings may reflect ethnic diversity (e.g., population 
stratification) or phenotypic heterogeneity. As a powerful tool for 
summarizing the results from different studies to estimate the 
major effect with enhanced precision, meta-analysis has generally 
been used in quantitative assessment of genetic variation and 
disease. Here we present the most comprehensive meta-analysis 
for the effects of E23K polymorphism oi KCNJll on T2D risk. 

Materials and Methods 

Literature search strategy and inclusion criteria 

To identify eligible literatures, we conducted a computer-based 
search of PubMed, Medline, EMBASE and ISI Web of Science 
databases without language restrictions. Studies pubhshed before 
the end of Mar. 2013 on T2D and the E23K polymorphism in the 
KCNJl 1 gene were retrieved. Search keywords combinations were 
"potassium inwardly-rectifying-channel, subfamily-J, member 11", 
"KCNJll", "Kir6.2", "type 2 diabetes", "type 2 diabetes 
meUitus", "T2D", "T2DM", "non-insulin-dependent diabetes 
meUitus", "NIDDM", "polymorphism" or "variation". The titles 
and abstracts were read to determine their relevance, and 
potentially relevant studies were retained for further evaluation. 
For retrieved articles, the fuU texts were carefully read to 
determine whether they meet the purpose of the present meta- 
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analysis. Furtlicnnore, the references of these studies were checked 
to identify other relevant publications. 

Eligible studies should meet following criteria: (1) focusing on 
the assoriation of the KCNJll E23K polymorphism with T2D risk 
(2) being case-control or cohort studies (3) diagnosis of T2D 
patient was confirmed j)atiiologically and (4) providing sufficient 
data for calculation of odds ratio (OR) with its 95% confidence 
interval (95% CI) and P-value. The major reasons for exclusion 
were (1) case-only studies (2) overlapping data (3) review papers. If 
more than two studies reported the same sample, only the study 
providing more information or latest published was selected. 

Data extraction 

The following information was carefully extracted from all 
eligible publications: the first author, year of publication, country 
of origin, ethnicity of subjects, study design, sample size, sex 
distribution among cases and controls, mean age and body mass 
index (BMI) of cases and controls, source of control, Hardy- 
Weinberg equilibrium (HWE) status in controls, genotyping 
method, number of genotypes in cases and controls. Two authors 
independendy assessed the articles for compliance with the 
inclusion criteria, and disagreement was followed by discussion 
until consensus was reached. 

Statistical methods 

The association of the KCNJll E23K polymorphism with T2D 
was evaluated by calculating a pooled OR and 95% CI for allele 
contrast (K vs. E allele), heterozygous (KE vs. EE) and homozygote 
(KK vs. EE). Then, we examined the association between the 
polymorphisms and T2D risk using dominant and recessive 
genetic models. The standard Qjstatistic test was performed to 
evaluate whether the variation between studies was due to 
heterogeneity or due to chance [7]. ORs were calculated 
according to the method of DerSimonian and Laird, and 95% 
CI was constructed by Woolfs method [8,9]. In addition, 
subgroup analysis was used to investigate potential sources of 
heterogeneity by stratified meta-analyses based on ethnic group, 
sample size (No. cases >1000 or <1000) and mean BMI of cases 
(<25, 25~30, or >30). Ethnic group was defined as East Asians, 
Caucasians (e.g., people of European origin), Indians and others 
(e.g., African American, Jews, and Arabian). Subsequendy, 
ethnicity, sample size, BMI, age and sex were analyzed as 
covariates in meta-regression to further investigate potential 
sources of heterogeneity. For family-based association studies, 
the transmission disequilibrium test (TDT) was used to analyze 
effect size of the polymorphism. In general, the OR was calculated 
from the ratio of transmitted alleles to non-transmitted alleles from 
heterozygous parents to affected offspring [10,11]. Combined 
effect size from both case-control and family-based association 
studies were calculated according to the method described 
previously by Lohmueller et al [12]. The Z-test was used to 
determine the significance of overall OR. 

We calculated the sample size required for 80% power with the 
summary OR estimated from each ethnic populations, assuming 
an equal number of cases and controls, risk allele frequency (RAF) 
in controls estimated from different ethnicity. Furthermore, 
population attributable risk (PAR) was calculated to get a 
comprehensive view of the impact of the E23K variant on T2D 
at population level. PAR was calculated by the following formula: 
(OR-l)/OR * risk allele frequency [13]. 

Egger's test and funnel plots were used to assess small studies 
effects [14]. Sensitivity analysis was performed by excluding one 
study at a time to assess the stability of the results. The type I error 
rate was set at 0.05 for two-sided analysis. AU of the calculations 



were performed using the ST ATA 10.0 (STATA Corporation, 
CoUege Station, TX) and SAS (version 9.1; SAS Institute, Cary, 
NC). 

Results 

Characteristics of included studies 

The literature search yielded 159 studies using keywords listed 
above. Figure SI shows the literature search and selection process 
for eligible studies (Figure SI). Finally, a total of 48 studies 
including 56,349 cases, 81,800 controls and 483 family trios, were 
retrieved based on the search criteria for T2D susceptibility related 
to the KCNJll E23K polymorphism [15-62]. In addition, almost 
all studies indicated that the distribution of genotype frequencies 
among the control groups were consistent with HWE. The 
detailed characteristics of all the included studies of this meta- 
analysis were summarized in Table 1. 

Meta-analysis results 

Overall, significant associations between KCNJll E23K poly- 
morphism and T2D were detected when all the eligible studies 
were pooled into the meta-analysis (Table 2). The overall result 
showed that the 23K allele of rs5219 polymorphism was 
significandy associated with elevated T2D risk with per-allele 
OR of 1.12 (95% CI: 1.08-1.17, /'<10"^ Figure 1). Significant 
increased T2D risks were also detected for heterozygous 
(OR= 1.09, 95% CI: 1.03-1.14, P<10"^) and homozygous 
(OR= 1.26, 95% CI: 1.17-1.35, P<10"^) when compared with 
wild type homozygous. Similar results stiU maintained using 
dominant and recessive genetic models (Table SI). When studies 
were stratified for ethnic populations, significant associations were 
also observed among East Asian and Caucasian populations with 
per-allele OR of 1.13 (95% CI: 1.08-1.17, P<10"-') and of 1.12 
(95% CI: 1.08-1.16, P<10"^) respectively. Significantly increased 
risks were also found for heterozygous and homozygous (Table 2). 
However, no such association was detected in Indian and other 
ethnic populations in all genetic models. In the subgroup analysis 
by sample size, significant associations were also observed for both 
large and small studies in all genetic models. When stratified by 
mean BMI of cases, statistically significant results were also 
observed for T2D cases with different BMI (Table 2). For two 
family-based association studies induding a total of 483 family 
trios, we failed to detect statistically significant evidence for the risk 
23K allele over-transmission from heterozygous parents to their 
T2D offspring (pooled ORtdt = 0.87, 95% CI: 0.72-1.05; 
P=0.14). 

Significant heterogeneity was found among the 46 included 
studies (P<10 *). Hence, meta-regression was further conducted 
to investigate the source of heterogeneity. In meta-regression 
analysis, ethnicity (P=0.79), sample size (P=0.61), mean age 
(P=0.36) of cases and controls (P=0.61), gender distribution in 
cases (P= 0.96) and controls (P= 0.30) did not explain a large part 
of the heterogeneity among the individual study. By contrast, 
mean BMI (P=0.03) explained about 11% of the heterogeneity. 

The 23K allele frequency of the rs5219 polymorphism varies in 
the control groups across different ethnic populations, ranging 
from 0.07 to 0.61 (Figure 2). In Caucasian controls, the K allele 
frequency was 0.40 (95% CI: 0.37-0.42), which was higher than 
that of East Asian controls (0.36; 95% CI: 0.34-0.38), Indian 
controls (0.34; 95% CI: 0.27-0.41). 2500 and 2200 case-control 
pairs will be required for 80% power to detect the risk allele 
among Caucasian and East Asian population respectively. The 
population attributable risk (PAR) of T2D related to E23K 
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Table 2. Results of meta-analysis for KCNJ11 E23K polymorphism and T2D risk. 



Sub-group 
analysis 


No. of 
studies 


K allele 










Heterozygous 






Homozygous 












OR (95°/oCI) 


flZ) 








OR (95%CI) 




m)' 


flO)" OR (95%CI) 




r\Q)' 




Ethnicity 












0.10 








0.06 






0.01 


Caucasians 


22 


1.12 (1.08-1.16) 


<10" 


5 


0.001 




1.09 (1.06-1.13) 


<10"^ 


0.13 


1.33 (1.18-1.50) 


<10"= 


0.0008 




East Asians 


14 


1.13 (1.08-1.17) 


<10" 


5 


0.02 




1.13 (1.05-1.22) 


0.0009 


0.02 


1.30 (1.18-1.42) 


<10"^ 


0.06 




Indians 


5 


1.06 (0.87-1.29) 


0.56 




<10"' 




1.00 (0.82-1.23) 


0.98 


0.01 


0.90 (0.69-1.18) 


0.45 


0.02 




Others 


7 


1.09 (0.97-1.23) 


0.15 




0.01 




0.98 (0.76-1.28) 


0.91 


<10"'' 


1.19 (1.00-1.43) 


0.05 


0.31 




Sample size 












0.32 








0.14 






0.04 


Large 


22 


1.12 (1.08-1.17) 


<10" 


5 


<10"^ 




1.12 (1.10-1.15) 


<10"'' 


0.12 


1.17 (1.07-1.28) 


0.007 


0.02 




Small 


26 


1.13 (1.07-1.18) 


<10" 


5 


<10"' 




1.10 (1.02-1.19) 


0.009 


0.002 


1.32 (1.20-1.46) 


<10"= 


0.0007 




Mean BMI of cases 












0.37 








0.20 






0.03 


<25 


12 


1.15 (1.11-1.21) 


<10" 


5 


0.15 




1.15 (1.08-1.22) 


<10"= 


0.21 


1.40 (1.28-1.55) 


<10"= 


0.74 




25-30 


25 


1.12 (1.06-1.19) 


<10" 


4 


<10"'' 




1.13 (1.06-1.21) 


<10"'' 


0.004 


1.18 (1.03-1.34) 


0.008 


<10"'' 




>30 


6 


1.10 (1.03-1.18) 


0.008 




0.07 




1.11 (1.04-1.19) 


0.001 


0.09 


1.16 (1.01-1.33) 


0.02 


0.10 




Total 


48 


1.12 (1.09-1.16) 


<10" 


5 


<10"^ 




1.09 (1.03-1.14) 


<10"= 


0.0001 


1.26 (1.17-1.35) 


<10"^ 


<10"= 





P(Z): Z test used to determine the significance of the overall OR. 

P[QlT: Cochran's chi-square Q statistic test used to assess the heterogeneity in subgroups. 
P(Q)'': Cochran's chi-square Q statistic test used to assess the heterogeneity between subgroups. 
doi:1 0.1 371 /journal.pone.0093961 .t002 



polymorphism was 4.6% overall, 4.4% for Caucasians and 4.5% 
for East Asians. 

Sensitivity analyses and publication bias 

The results of sensitivity analysis confirmed the significant 
associations of the KCNJll E23K polymorphism with T2D risk, 
and no single study influenced the overall OR qualitatively (Figure 
S2). The Egger's test and funnel plots indicated no publication bias 
for the association of KCNJll E23K polymorphism and T2D 
(Figure S3; Egger test, /^0.05). 

Discussion 

Limited statistical power of relative small sample size is a 
common problem in genetic association for individual T2D 
studies. Therefore, sufficient sample power is necessary in 
deciphering genetic architecture of T2D, but it is sometimes very 
difficult for a single study to collect enough amounts of data to 
reach a reliable conclusion. By pooling of data from individual 
association studies, meta-analysis is an effective approach of 
increasing the sample size under investigation, thus enhancing the 
statistical power for the estimation of genetic effects. 

Our results indicated that the rs5219 polymorphism oi KCNJll 
is a risk factor for developing T2D. In the subgroup analyses by 
ethnicity, we found the rs5219 polymorphism was associated with 
T2D among East Asians and Caucasians, but not Indians or other 
ethnic populations. Of note, different ethnic populations were 
pooled in the other ethnic group and only a few studies were 
available in the subgroup, so the result must be interpreted with 
caution. There are several other possible reasons which may 
account for such dififerences. First, T2D is a complex disease and 
different genetic backgrounds may cause the discrepancy since the 
distributions of the risk- association alleles in KCNJll were 
different between various ethnicities. The K allele frequency of 
the rs5219 polymorphism was ~36%, ~40%, and ~34%, among 
East Asians, Indians and Caucasians populations, respectively. 
Such result could also be due different Knkage disequilibrium (LD) 



pattern of the polymorphism and nearby causal variant among 
different ethnic populations. Moreover, inter-individual difference 
like age, sex, dietary intake of nutrients, in addition to phenotype 
heterogeneity, such as years from onset and severity of the disease 
may also explain the discrepancy. Furthermore, study design and/ 
or small sample size or some environmental factors may also affect 
the results. Therefore, more studies are needed to further validate 
the effect of the polymorphism on T2D risk among difference 
ethnic populations. 

The KCNJll gene has attracted considerable attention as a 
promising candidate for T2D based on its position and its function 
as a key factor in the regulation of glucose-induced insulin 
secretion, since normoglycemic lysine carriers are shown to 
consistently display a defect in insulin secretion [21,63,64]. 
Functional studies suggested that the KK genotype might induce 
a critical inhibition of glucose-induced insulin release from 
pancreatic fi-cells [65]. Furthermore, the KCNJll E23K variant 
was found to be associated with glucose intolerance and 
conversion from impaired glucose tolerance to T2D among 
Caucasians [66,67]. Previous studies indicated that the E23K 
variant is functional by affecting in vitro properties of KATP 
channel via increasing the threshold of ATP concentration for 
insulin secretion [65,68]. 

The distribution of the E23K variant in controls across various 
studies showed global variation (Figure 2), suggesting the 
possibility of population stratification. However, empirical evi- 
dence indicated that well-designed population-based association 
studies can keep the effects of population stratification to a 
minimum [69]. Almost all the population-based genetic associa- 
tion studies included in the present meta-analysis were well- 
designed by recruiting cases and controls from the same 
geographic region and ethnicity, which may help to reduce the 
effects of population stratification. Moreover, the effects of 
potential population stratification in any individual study may be 
in a random direction, so that one individual study with a small 
amount of stratification should have very limited effect on the 
overall results [70]. By combining results from TDT studies, which 
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Figure 1. Forest plot from the meta-analysis of T2D risl< and KCNJ11 rs5219 polymorphism using random-effects model. 

doi:1 0.1 371/joumal.pone.0093961 .gOOl 
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Risk allele frequence stratified by ethnicity 



Frequence 
1.0 




East Asians 



Caucasians 



Indians 



Others 



Figure 2. Frequencies of the 23K allele of KCNJ11 rs5219 polymorphism among controls stratified by ethnicity. The "□" represent 
outlier. 

doi:10.1371/journal.pone.0093961.g002 



are robust to potential population stratification, we failed to detect 
an over-transmission of the 23K allele from heterozygous parents 
to their T2D offspring. Given the small number of studies and 
relative small sample size, the combined results from TDT studies 
should be interpreted with caution. 

Power analysis revealed a hint in design for future association 
studies. So far, most association studies have included at the most 
several hundred subjects and results from these studies should be 
treated with caution as for limited statistic power to reach a 
reliable conclusion. According to our power analysis, powerful 
association studies on the E23K polymorphism and T2D risk may 
need several thousand individuals. Compared with candidate gene 
approach, genome wide association study (GWAS) with large 
sample size and unbiased to genomic structure is a powerful 
approach in susceptible gene identification for T2D. Recently, 
Tsai et al. [71] reported a two-stage genome-wide association 
(GWA) conducted in Chinese and identified KCNJll as a risk 
region in T2D susceptibility which was in line with the results the 
present meta-analysis. 

In comparison with the previous published meta-analysis 
[63,66,72-74], the current study included more than ten times 
as many cases as the earlier meta-analysis. In addition, we assessed 
the effect of risk allele and T2D using various genetic models and 
reached consistent results. Furthermore, we systematically ex- 
plored potential sources of heterogeneity across studies and the 
possibility of publication bias. In addition to combine those newly 
published data, the present study also statistically joined popula- 
tion-based and family-based genetic association studies into a 
single meta-analysis, which allowed us to enhance the power of the 
meta-analyses and also establish a comprehensive picture of the 
relationship between KCNJll and genetic susceptibihty to T2D. 

Heterogeneity among pooled studies is a frequently encountered 
issue in meta-analyses. Of note, our meta-analyses joined 
population bases association studies and family based association 



studies, which could enhance the heterogeneity. Hence, we 
performed subgroup analysis and meta-regression to identify 
potential source of between-study heterogeneity. However, they 
revealed that only BMI of controls could explain a small part of 
significant heterogeneity between studies. There may be a number 
of possible underlying reasons. Firstly, results from case-control 
studies may differ because of ethnic diversity (e.g., variation in 
allele frequencies) and geographic variation. Secondly, variations 
in methods of sample ascertainment and diagnosis may also 
contribute to such inconsistence. Thirdly, environmental factors, 
such as alcohol drinking, smoking behavior, obesity might also 
result in variability between studies. 

To conclude, this may be the most comprehensive meta-analysis 
of KCNJll and T2D. Our results suggest a modest but statistically 
significant effect of the 23K allele of rs5219 in susceptibility to 
T2D, particularly in East Asians and Caucasians. More work wiU 
be required for future association studies, especially those which 
are properly powered, effectively control for confounding factors, 
employing family-based design. Moreover, gene-environment and 
gene-gene interactions should also be taking into consideration for 
future studies. 

Supporting Information 

Checklist SI 

(DOC) 

Figure SI Flow chart of literature search for studies 
examining KCNJll rs3219 polymorphism and risk of 
T2D. 

(TIF) 

Figure S2 Result of sensitivity analyses. 

(TIF) 



PLOS ONE I www.plosone.org 



7 



April 2014 I Volume 9 | Issue 4 | e93961 



Type 2 Diabetes Genetics 



Figure S3 Begg's funnel plot of rs5219 polymorphism 
and T2D risk. 

(TIF) 

Table SI Results of meta-analysis for KCNJll E23K 
polymorphism and T2D risk using dominant and 
recessive genetic models. 

(DOCX) 
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